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We discuss the excitation spectra around the Dirac node on a surface of a three-dimensional 
topological insulator. By using the diagrammatic expansion, we show that the coupling of an 
electron with the gauge field in the presence of impurity scatterings opens a gap around the Dirac 
node. The results are consistent with a recent experimental finding by T. Sato, et al. [Nature 
Phys. 7, 840 (2011)]. We also discuss the consistency between the present results and the bulk- 
boundary correspondence of the topological insulator in the presence of time reversal symmetry. 
The conclusion can be applied any two dimensional massless fermions. 
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I. INTRODUCTION 

Properties of metallic states at the surface of topolog- 
ically non-trivial materials are a hot issue in the con- 
densed matter physics. A three-dimensional(3D) topo- 
logical insulator (TI) hosts a metallic state with linear 
dispersion at its surface as a result of the topologically 
non-trivial nature of the wave function in bulk insulat- 
ing region^— . According to the shape of the dispersion, 
such metallic state is called Dirac cone. A point in the 
momentum space at which the upper cone and the lower 
one touch with each other is called Dirac node. The 
bulk-boundary correspondence suggests the presence of 
a single Dirac node at the surface, which reflects a topo- 
logical number Z2 = 1 in the bulk insulating region. The 
Dirac node is considered to be fragile against perturba- 
tions which break the time reversal symmetry (TRS) be- 
cause Zi is well defined only in the presence of TRS. In 
fact, attaching a ferromagnetic insulator onto a TI^£ or 
introducing magnetic impurities into a Td£ remove the 
Dirac node from the surface state. On the other hand, 
the Dirac node is believed to be robust against perturba- 
tions which preserve TRS. The Dirac node remains even 
in the presence of normal impurity scatterings, which is 
an important feature of the surface state from theoretical, 
experimental and applicational view points. 

A recent experiment^, however, has reported the 
gapped excitation spectra in the single Dirac cone at the 
surface of TlBi 2 (Se x Si_ x )3. Starting with a topological 
insulator TlBi2Se3, the gradual substitution of Selenium 
by Sulfur enables to make a series of insulators ended 
with a topologically trivial insulator TlBi2S3. The insu- 
lators should be topological and should have a gaplcss 
single Dirac cone at their surface in the doping range of 
0.5 < x < 1.0. The experimental results of the angle re- 
solved photo emission, however, clearly shows the gapped 
excitation spectra in the corresponding doping range. At 
present, we do not have any reasonable argument which 
explains the discrepancy between the theoretical predic- 
tion and the experimental results. 

Motivated by the experiment^, we theoretically try to 



make clear a mechanism which generates the gap around 
the Dirac node in the presence of the TRS. We con- 
sider the massless fermion in two-dimension which cou- 
ples with the gauge field in the presence of normal im- 
purity scatterings. The self-energy due to the impurity 
scatterings and the coupling with the gauge field is calcu- 
lated within the lowest order of perturbation expansion. 
The results show that the real part of the self-energy for 
the upper Dirac cone has the opposite sign to that for 
the lower cone. As a result, two-dimensional fermion has 
gapped energy spectra around the Dirac node. We con- 
clude that the interplay between the impurity scatterings 
and the dynamical electron-electron interactions via the 
gauge field generates the gap. The results suggest the 
absence of gapless surface state even in the presence of 
TRS. We will show that our conclusion is not inconsistent 
with the bulk-boundary correspondence of a TI. We also 
conclude that the translational symmetry is necessary to 
preserve the gapless surface state in real TIs. 

This paper is organized as follows. In Sec. II, we ex- 
plain our theoretical model. The self-energy due to the 
coupling with the gauge field in the presence of the im- 
purity scatterings is calculated in Sec. III. We discuss a 
relation between theoretical results and an experimental 
one in Sec. IV. The conclusion is given in Sec. V. 

II. THEORETICAL MODEL 

Let us consider the Lagrangian density which describes 
the two dimensional massless Dirac fermions coupling 
with the gauge field, 

E ^fe-^'W-J^ (!) 
3=0,1,2 

where er° is the 2x2 unit matrix, er 3 for j = 1 — 3 
are the Pauli matrices, i/> and A are the field operators 
of a surface massless fermion and a U(l) gauge boson. 
In Eq. ([1]), Po = Po represents energy measured from 
the Dirac node and pj = Pj/v for j = 1 — 2 represents 
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the momenta multiplied by the speed of light c where 
v = vp/c is a small constant describing dimensionless 
velocity with the Fermi velocity vf- Throughout this 
paper, we use the units of h = c = 1. The field tensor 
F^ v is denned by 



(2) 



with <9 M = d/dx^. We consider the Dirac propagator 
defined by 



Go(p) = 



Po<?" 



vp ■ cr 



Po- 



tt 2 p 2 



■ is 



(3) 



where we take a short- hand notation p = (po,p). Eq. ((3]) 
represents the particle propagator for positive energy 
Pa > and the hole propagator for negative energy 
Po < 0. We employ the Feynman gauge in which the 
propagator converges at the spacial infinity Xj — > oo with 
j = 1 — 2. The propagator of the gauge field at energy 
qo and momentum q is given by 



qo - q + i£ 
gpv =diag{l, -1, -1}, 



(4) 
(5) 



where fi and v are 0, 1, and 2. As shown in Eq. ([1}, 
the vertex of the coupling between the electron and the 
gauge field is described by — ea^ in Heavisidc units. We 
assume that the impurity potential is spin-independent 
and is represented by the delta-function, 



U i (r)=^2va°6{r-r j ), 



(6) 



where rj is the position of an impurity and u is the 
strength of a single impurity potential. We assume a 
constant number density of impurities n,. 



(a) 




(b) 




FIG. 1. The Feynman diagrams under consideration. Figure 
(a) represents the fermion self-energy coupled to the gauge 
boson. The wavy line denotes the propagator of the gauge 
boson. Figure (b) is Dirac fermion propagator in the presence 
of impurity scatterings within the Born approximation. 



The details of calculation are shown in Appendix [B] 
where we estimate the self-energy in Fig. [TJa) . Namely 
we calculate T,^(p) in Eq. (jS]) with using the full Green 
function G(p—q) in Eq. (|B1|) instead of the bare Go(p—q) 
in Eq. ([3]). In the limit of small momenta, the real part 
of the self-energy has an asymptotic form as 



Re 



E<»>(p) 



Po& 



p ■ cr 



(9) 



We note that the result does not contain 7. Therefore the 
impurity scatterings add only negligible corrections to 
Eq. @ as shown in Appendix [5] From pole of the scalar 
denominator [Go(p) _1 — this self energy cannot 

open gap for massless Dirac fermion, 



III. SELF-ENERGY 

At first, we estimate the self-energy by the impu- 
rity scatterings alone. The impurity potential mixes the 
states with different momenta. It, however, does not mix 
the states with different energies. Therefore impurity 
scatterings itself cannot open the gap at the Dirac node. 
The self-energy within the Born approximation is repre- 
sented by 

E «) = (7) 

where 7 = iru 2 ni/v 2 > is the dimensionless expan- 
sion parameter smaller than unity (See also Appendix 
IA")) . The result has the general form of the massless Dirac 
fermion^. 

Secondly, we calculate the self-energy of a Dirac 
fermion coupled with the gauge field which is given by, 

= ££(-e^)G (p-g)(-eOJV(9)- ( 8 ) 

q fiv 



Po=± (v\p\ + ire 2 ) - 7re 2 sgn(p ) 

= ±v\p\. (10) 

The two terms stemming from the self-energy cancel each 
other out, which reflects the covariance of the gauge field. 
As a consequence, the coupling with the gauge field does 
not change the energy spectra of massless Dirac fermion. 
Neither the coupling with gauge field alone nor the im- 
purity scatterings alone open the gap at the Dirac node. 
The former changes the energy of a fermion but preserves 
the translational symmetry. On the other hand, the lat- 
ter breaks the translational symmetry and changes the 
momenta of a fermion. But it preserves the energy of 
fermion. To have a gap at the Dirac node, we need scat- 
tering processes which change the momenta and the en- 
ergy of a fermion at the same time. 

Finally, we consider a self-energy which represents the 
interplay between the coupling with the gauge field and 
the impurity scatterings. The scattering processes are 
diagranimatically described in Fig.[^a). Such self-energy 
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FIG. 2. The Feynman diagrams we consider in the text. 
Fig. (a) is the fermion self-energy consisting of the gauge boson 
and the impurity scatterings. Fig.(b) represents the vertex 
function including the impurity scatterings. 

is given by 

E(p) = J2 ^{p' ,q)G{p-q){-ea v )D^ u {q), 

(11) 

^(Po,q) = 7^ G(p)(-«T")G(p - 0), (12) 
p 

where r M (po, q) is the vertex function shown in Fig.[2Jb). 
The vertex function can be derived from the Ward- 
Takahashi's identit y 11 ' 12 with the self-energy of impurity 
scatterings (Eq. J7])), 

(p> - p)^ + TMp')) =P'- -{p- E(«)(p)) . 

(13) 

The vertex function results in 

r"(p,p') = -ie-fS ", (14) 

where 8^ v is the Kronecker's delta. The correction at 
another vertex of gauge field (Fig. [5] (a)) give the same 
contribution. Thus the self-energy in Eq. (JTT|) becomes 

S(p) = 2z 7 e 2 ^G(p- g ) 2 1 (15) 

The summation of q can be done in the same method 
for deriving Eq. ([8]). By introducing a parameter xq = 
v 2 p 2 /po 2 , we obtain the expression of the self-energy in 
the limit of xq — > as 

lim Efo.p) =2^e 2 7 2 ((sgn(pr>° + p ■ a) , (16) 

where p is the unit vector in the direction of p. This 
is the central result of this paper. The result does not 
depend on either the magnitude of the energy or that of 
the momentum. From the pole of the scalar denominator 
of [Go(p) _1 — S(p)], we obtain the energy dispersion of 
perturbed fermions at |p| — > as 

Po = ± {v\p\ + 2?re 2 7 2 ) + 27re 2 7 2 sgn(p ) 

= ± (v\p\ +47re 2 7 2 ). (17) 



The result shows the gapped energy spectra around the 
Dirac node. We note that both the impurity scatterings 
and the coupling to the gauge field preserve the TRS. 



IV. DISCUSSION 

According to the features of the spectra in Eq. (|17l) . 
the magnitude of the gap should be independent of tem- 
perature. Our results are consistent with the gapped 
spectra found in the experiment by T. Sato et al£. An- 
other experiment^ 3 - has reported the gapless spectra in 
the similar situation. The energy resolution of experi- 
ment in Ref. seems to be better than that in Ref. [l3|. 

We discuss effects of electron-electron interactions via 
the gauge field on the excitation spectra around the Dirac 
node. Because of the dynamics of the gauge field, the in- 
teraction potential depends on frequency As a result, 
the interplay between the impurity scatterings and the 
interaction via the gauge field opens the gap at the Dirac 
node. On the other hand, in condensed matter physics, 
the electron-electron interactions are usually taken into 
account within the static approximation. Namely, the 
Coulomb interaction acts two fermions at an equal time. 
Therefore such static interaction potential does not de- 
pend on the frequency. The two different types of inter- 
actions in the presence of impurity scatterings may affect 
the excitation spectra in two different ways. In one- 
dimension, sufficiently strong static interactions cause 
the gap^. In two-dimension, however, the excitation 
spectra have been believed to remain gapless even in the 
strong static interactions. It is necessary to consider the 
electron-electron interactions via the dynamical boson 
field to have the gap at the Dirac node in two-dimension. 
This argument is also important when we consider topo- 
logical protection of the gapless surface state by the bulk- 
boundary correspondence, which we discuss next. 

Finally, it is necessary to discuss a relationship be- 
tween our results and the bulk-boundary correspondence 
of three-dimensional TI in the presence of the TRS. The 
Zi invariant is well defined in (1+3) dimensional fermion 
space of Tli. The bulk-boundary correspondence works 
well to understand the existence of the two-dimensional 
gapless state on its surface. In this paper, however, we 
consider a situation in which the fermion in TI naturally 
couples with the electromagnetic field. In such situation, 
the whole physical space consists of the fermion space and 
the gauge boson one. Namely, the (1+3) dimensional 
fermion space is no longer closed independent physical 
space. At present, it is unclear if it is possible to apply 
the topological classification defined solely in the fermion 
space to the integrated physical space of the fermion and 
the gauge boson. Even in the integrated space, it may 
be still possible to consider the topological characteriza- 
tion in the fermion space only. In such case, the theory 
would be possible to define some topological invariants in 
a similar way as the weak topological invariants^. There- 
fore the predictions by the topological theory^ " 3 ' 16 ' 17 in 
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(1+3) dimensional fermion space are not always valid. 
This makes the background of the gapped energy spec- 
tra appearing at the topologically protected surface state. 
The coupling of fermion with another degree of freedom 
sometimes changes a topologically non-trivial phase to 
a topologically trivial one. Our result demonstrates an 
example of the story. 



V. CONCLUSION 

In conclusion, we have studied effects of the gauge 
field and the impurity scatterings on the energy spec- 
tra of massless Dirac fermion at the surface of a three- 
dimensional topological insulator. By using the per- 
turbation expansion, we have calculated the self-energy 
around the Dirac node. The interplay between the im- 
purity scatterings and the interaction with gauge field 
leads to the gapped energy dispersion around the Dirac 
node. Our results suggest that the translational invari- 
ance is also necessary for the surface state to be gapless 
in three-dimensional topological insulators. 
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Appendix A: Self-energy due to impurity scatterings 

We derive the self-energy of impurity scattering alone. 
Within the Born approximation, the self-energy is given 
by, 



J po — V z p z + IE 



(Al) 



The term proportional to momentum p vanishes after 
integrating of p. The self-energy is calculated as 



L v ' =U 7 



d 2 p 



Po& 



=u riipaa 



(v2tt) 2 Po - i ' P + ie 

.0 



2 I 2 

7 =iru rii/v 



d 2 p f 00 da 

WW L H) 

o 



cxp[ia(p — v p + ie)] 

(A2) 
(A3) 
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The result is pure imaginary and gives the inverse of the 
life time in agreement with a previous result^. We intro- 
duce a dimensionless constant 7 which contains material 
parameters. 



Appendix B: Self-energy due to coupling with gauge field 

Let us calculate the self-energy due to coupling with the gauge field represented by Fig. HJa). The propagator for 
a Dirac fermion in the presence of impurity scattering is represented by 

(1 + ij) po — up 
By using this propagator, the self-energy in Fig. (Ha) is given by 

^^(-e/lGfp-?)^)^)?) (B2) 

q )XV 

The self-energy is expressed in the parametric integral representation^ as, 

/oo j j2_ roc poo 1 
7^=4/ / daida 3 (-e<^)-((T d Xo -a-a a ,)(-ea v )exp[iS], (B4) 
-00 (V27r) 3 Jo Jo * 

S =ai((l + i'jfipo - qof - v F 2 {p - qf + ie) + a 2 (q 2 - q 2 + ie) + x Q (l + 17) (po - qa) - vx ■ (p - q), (B5) 
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where x M and ctj are dummy variables. The integral with respect to the momenta can be done by the Fresnel integral. 
The results become 



e 2 lim 



oo poo 



da\da 2 



x cxp 



*->°Jo Jo {a x v 2 + a 2 )y / i(a 1 (l + ij) 2 + a 2 ) * 

a\a 2 {\ + ij)po 2 aia 2 p 2 a 2 (l + i^)p x vx ■ p 



i < i(ai + a 2 )e + 



«i(l + 17) + a 2 a\V 2 + a 2 «i(l + ij) 2 + a 2 a^v 2 + a 2 



Next we introduce a integration variable p = ai + a 2 and change the integration variables as on — > pat. In this way, 
we obtain 



E^fp) 



da\da 2 5(1 — a\ — a 2 ) 



Jo a i y2 + a 2 + i~f) 2 + a 2 ) \ai(l + i-f) 2 + a 2 

a\a 2 (l + i-f) 2 



a 2 

2 ; vp ■ a 



aiV + a 2 



r d p 




/ ^=exp 


'->( 


/o V«P 





i(ai(l + i-f) 2 + a 2 ) ° i{a\v 2 + a 2 ) 
By applying the condition » 2 <Cl and introducing a variable xo = v 2 p 2 /po 2 , the self-energy becomes 

da 



a x a 2 ^ 



1+27 Vp ■ <T 

+ i7) 2 po 2 (l - a) - (1 + (2z7 + 7 2 )a)v 2 p 2 V 1 + 2z7a PoCr 1-Q 

1 



'Per 



p(l-*7) 



da 



7(21 + 7) 1-a 
1 + . x 



^((1 - x ) - a) V 2(1 + i 7 ) 2 1 - a - x J \ 1 + 2i 7 a 



1 + 17 n vp ■ a 
Po° ~ 

1 — Q' 



(B6) 



(B7) 
(B8) 



We neglect terms proportional to xq because we focus on the limit of xq ~ in the followings. The results are 



^ 9 \p) 



da 



Po 2 Wi-i iy/a(a - (1 - x )) 



1— XQ 



do 



\/a((l - Xq) - a) 



(1 + i-f — 2i'ya)po<J 



' o- 

1 - a 



(B9) 



VPo' 
(1 + ij)po<J° 
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We reach the self-energy as, 



£ (9) (p) 



(1 - 27) + i7)p o" lo{ 



+i7r7Pocr xo - i — log(l + x ) - 7r 



1 + x + 2^/x^ 

1 - x 
vp ■ cr^ 



+ irpoa - 27p o' v / ^o + 272?oc°(l - x ) log 



1 - y/xp 



1 + ^/x~o~ 
(Bll) 
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Using xo = v p /po <C 1, it becomes 



SW(P) c - e 2 (l - i 7 ) - ^ - 2 7^^4 + 2 7 (1 - xo) log (± 

V VPo VP V-Po \ x 



1 - V^o \ Po^ 



-z(l + 17J — log - h iit^xq- 



1 - x 



ilog(l + x ) 



p ■ a 

Vp 1 , 



(B12) 



Finally by putting 7 — > 0, it is possible to obtain the self-energy due to the gauge field alone. 
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